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A Graphical Method of Drawing Trajectories for High- Angle 

Fire, 

By W. E. Dalby, F.R.S. 
(Received January 25, 1917.) 

In May, 1915,* I described a graphical method for finding the time of 
flight, the range, the angle of elevation, in fact, the elements of a trajectory, 
having given the curve of resistance as a function of the velocity and 
assuming that the density of the air was constant at all points along the 
trajectory. The method is therefore only applicable to trajectories of low 
altitude. In this paper I describe a graphical method for obtaining a series 
of points on a trajectory of high altitude, taking into account the variation 
of the density of the air with height. 

Let A, fig. 1, be a point on a trajectory where the velocity is given and 
equal in direction to <^^ and in magnitude to V^, and whose position is given 
by the co-ordinates X and Y. The problem is to find, the co-ordinates ^, y, 



with reference to the point A, of a second point B, on the trajectory, taking 
into consideration the actual resistance of the air to the motion of the shell, 
it being assumed that the variation of the density of the air is known as the 
shell rises or falls in height ; and to find the angle cf)^ defining the direction 
of motion and the magnitude of the velocity at the point B in the trajectory. 
Providing that two curves are given, the problem may be solved approxi- 
mately by a process of step-by-step integration, using the graphical method 

^ * Roy. Soc. Proc.,' A, vol. 92, p. 329 (1916). 
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explained in my paper* to obtain time-velocity and time-distance curves. 
These two curves are : — 

(1) A curve showing^, the resistance in air of standard density of a shell 
of standard form fired under standard conditions, the resistance being shown 
as a function of the velocity (such a curve is shown in fig. 3). 

(2) A curve showing the tenuity r of the air as a function of the height 
(as shown in fig. 7, p. 344). The tenuity r at height h is the ratio between 
the density of the air at height /z. and the standard density at sea-level. 

According to the official text-books the resistance to motion through the 
air is then proportional to rj^, and for a shell of diameter d the resistance is 
equal to JcacPrp, where y^ is a coefficient to allow for the shape of the shell 
and o- is a coefficient to allow for the steadiness, a technical term used to 
include such variations as may be produced by irregularities in the actual 
firing. The resistance expressed in this way is also in accord with the 
method given in the official gunnery text-book. 

Sohition of Problem. 
The given velocity at the point A in the trajectory (fig. 1) has components 

Va sin <^A vertically, 
Va cos <^A horizontally. 

The resistance in the trajectory corresponding to the velocity V^ is found 
from the curve fig. 3 to be ^^. The corresponding component resistances at 
A are :— 

^A sin ^A vertically, 

Pa cos ^a horizontally. 

Then, putting E for the constant had?, the resistances against which the 
shell rises vertically are the air resistance Erp sin ^, and the weight of the 
shell w, giving a total resistance of (Erp sin ^ + '^o). When the shell is rising 
the two terms are of the same sign, but when falling they are of opposite 
sign. 

The resistance against which the shell moves horizontally is Erp cos ^. 
The equations of motion are : — 

~~(Etj9 sin (j) + w) =^!^ -- vertically, (1) 

Q do 

■— (Erjj COS <^) = - -- horizontally. (2) 

^ dt 

^- 'Eoy. Soc. Proc.,' A, vol. 92, p. 329 (April 1, 1916). 
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The formal solutions of these equations are 



t 



w 



\ 



B g (Erp sin <^ + w) 

dv. 



• dv, 



%v 



(3) 



(4) 



B g (Etp cos <^) 

In these expressions jp is a function of the velocity and r is a function of 
the height and is independent of the velocity. The equations cannot, 
therefore, be integrated as they stand, since t cannot be expressed as a 
function of v. The graphical method explained in my previous paper may, 
however, be applied step by step to get a solution in the following way > — 

(1) Assume a value of the vertical velocity at B near enough to the given 
vertical velocity, 'yAsinc^A, at A to make the influence of the change of 
tenuity on the resistance negligibly small. Let v^ stand for the assumed 
value of the vertical velocity at B. Plot a curve representing the quantity 



w 



_ (5) 

^(ETAjpsm^A + t^) 

between the values of the given vertical velocity at A and the assumed 
value at B. In plotting this curve take r constant and equal to the 
tenuity at the height A. The resistance ;p is the only variable, and its 
value at a few points along the curve and the final point corresponding to 




the assumed velocity is to be found by dividing the velocities selected by 
sin ^A in order to reduce them to actual velocities in the trajectory, and then 
the corresponding values of j? -are taken either from a Table, or from a curve 
like the curve in fig. 3, drawn to a large scale. The error in the resistance 

2 D 2 
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caused by assuming that ^a remains constant during the motion over this 
portion of the trajectory is negligibly small. 

(2) From the curve representing equation (5) above, construct the corre- 
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spending parts of the time-velocity and the time-height curves by the 
graphical method described in my previous paper. 

(3) From the diagram obtained in (2), measure off the height corre- 
sponding to the assumed vertical velocity at B, namely, v^, and also the 
time-interval, and call the time-interval T. Equation (3) has now been 

/•A 

solved, and also its derivative, h ™ I vdt, on the assumption of constant. 

Jb 

tenuity. The solution gives— 

(a) The vertical step in height y, 

(h) The time-interval T corresponding to y. 

We now go to equation (4), solving it in a similar way. The link 
connecting the two equations is the time-interval T. Elements of the 
time-velocity and the time-distance curves are constructed from assumed 
values of the horizontal velocities and an appropriate value of the tenuity,. 
and then points are found on these elements corresponding to the time- 
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interval T, fixing at once the corresponding values of the horizontal step x 
and the horizontal velocity at the end of the step. 

(4) Assume a few even values of the horizontal velocities, divide them by 
cos (jf)A to get the corresponding values of the velocities in the trajectory, 
from which values of p are found for use in the expression 

(6) 



g (Etip cos (jf)A.) 

The value of n to be used in this expression is the average value during 
the step y. That is n = (ta + tb)/2. Values of r^. and tb can be found by 
calculation, or from a curve similar to that shown in fig. 7, where the 
tenuity is plotted as a function of the height. Generally, two values of (6) 
will be sufficient to determine the element of the auxiliary curve. 

(5) From the curve found in (4), construct the corresponding parts of the 
time-velocity time-distance curves by the general graphical method described 
in my previous paper. 

(6) From the curves found in (5), measure off — 

(c) The distance corresponding to the time-interval T. 
{d) The velocity at the end of the time-interval T. 

Then (c) is the horizontal distance x and {d) is the horizontal velocity at B. 

(7) It remains to determine the direction of the tangent at B. Let the 
tangents to the trajectory through the points A and B (fig. 2) meet in E. 
Draw the chord AB. Then the angle EAB is equal to ^a — ^- And it may 
be assumed, with negligible error, that the angle EBA is also equal to this. 
Then the exterior angle i is equal to 2 ((J^a — O), and the angle <j[)b, which the 
tangent at B makes with the horizontal, is given by 

(^B = </>A~2((^A~^). (7) 

The angle in this expression is to be calculated from 

tan = y/x, 

using the values of y and x found respectively in (a) and (c). 

(8) The velocity in the trajectory at B is given alternatively by 

the assumed vertical velocity at B 



V 



B 



or V 



B 



sin ^B 
the determined horizontal velocity at B 



cos ^B 

These two should substantially agree and together furnish a check on the 
accuracy of the work. 
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(9) The initial conditions at B are now : — 

Direction of the tangent <^b- 
Velocity along the tangent Vg. 
Tenuity tb and the air resistance j)B' 

All the data are therefore known for a fresh application of the method to 
find a point C on the trajectory. 

Eepeated applications of the method from point to point carry the 
trajectory up to the vertex, and then down again to any level which may 
be assigned. About eight points are sufficient to carry up to the vertex of 
the trajectory, and about six points bring it down to the ground again. 
Height to the vertex, and the range on the level or on any inclination, 
together with the time of flight, are then determined. 

I made a first application of the method to the case of a shell weighing 
380 lb., fired from a 9*2-inch gun elevated to 40°, taking a muzzle 
velocity 2375 feet per second, the data corresponding with one of the 
rounds fired from the 9*2-inch wire-wound gun used in the "Jubilee Trials'' 
on June 12, 1888. The results obtained were : range 23,000 yards, and 
time of flight 70*65 seconds. The actual range and time of flight recorded 
in the reports of the trial are: range 20,236 yards, time of flight 
63*5 seconds. 

A study of the original report shows that there were difficulties 
which probably reduced the muzzle velocity. Allowing for this probable 
discrepancy between the muzzle velocity assumed, namely, 2375 feet per 
second, and the probably lower value which was actually obtained on the 
range, the agreement is good. 

It may be noticed that, if the calculated range is reduced in the ratio of 
the observed time of flight to the calculated time of flight, the result gives a 
range very nearly equal to the observed range. Thus 

23000x63-5 on^^o i 

___ — --— -— = 20,672 yards. 

70-65 "^ 

The trajectory defined by 13 points found in the way described above, is 
shown in fig. 8 (p. 345). From this it will be seen that the greatest height is 
6215 yards, and the striking angle or angle of descent is 46^°. 

The method of drawing a trajectory is illustrated in the following example 
by giving in detail the work of determining a point on this particular 

trajectory. 

Uxample. 

A shell weighing 380 lb. is fired at an angle of elevation of 40^ with 
the horizontal, from a gun 9*2 inch diameter. The muzzle velocity is 
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2375 feet per second. The product of the coefficient of shape, h, and the 
coefficient of steadiness, a, is 0*9. The vahie of the coefficient E is therefore 
0-9 X 9-22 = 76-2. 

Plot the trajectory, and from it find the height of the vertex above the 
level, the time of flight, and the range. 

The resistance /j is given by the curve fig. 3. The tenuity t is given as a 
function of the height by the curve fig. 7, and is also indicated in fig. 5. 
The sections following are numbered to correspond with those included in 
the explanation of the method given above. 

The point A, fig. 4, is, for the first interval of the trajectory, the origin of 
the trajectory. The value of the tenuity r is unity, and the value of ^ 
corresponding to the muzzle velocity of 2375 feet per second is 21*03 lb. 

At A the vertical velocity is Va sin c^a — 2375 sin 40° = 1525 feet per 
second. 

At A the horizontal velocity is Va cos ^a = 2375 cos 40° = 1820 feet 
per second. 

The Vertical Component y, — (1) Assume a vertical velocity at B of 
1400 feet per second. 

The value of the quantity (5) is calculated from 

380 11-8 



32-2 (76-2 X 1 X p sin 40° + 380) 76-2 xp sin 40°+ 380* 



(8) 



At A, where the velocity is 2375 feet per second, p = 21*03. 

At B, where the assumed vertical velocity is 1400 feet per second, the 
corresponding velocity in the trajectory is 1400/sin 40° = 2180 feet per 
second, and the corresponding value of p is 18*49. 

(2) The values of the quantity (5) at A and B respectively are 0*008525 
and 0*009345. These are set up at the corresponding vertical velocities at 
Or and 1400—/ (fig. 5). It is unnecessary to calculate intermediate values 
for the interval, so that r and / are joined by a straight line. 

The area (a) represents the time taken for the vertical velocity to fall 
from the initial value, 1525, to the assumed value, 1400, and is, in fact, the 
solution of the differential equation (3) between the limits of the vertical 
velocities at A and B. The distance 1400— m represents this area, and m 
is therefore a point on the time-velocity curve, and Om is the first element 
of the curve. 

Again the area (b) represents the vertical distance moved through in the 

time 0^. The length tn is set out to scale to represent this distance, so that 
^ is a point on the time-height curve, and O^^. is the first element of the 
curve. 
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A vertical through n cuts the tenuity curve in q at the point where its 
value has fallen, to 0*956. 




(3) From this diagram (fig. 5), therefore, 
0?^ = 1'116 seconds is the time interval T, 
Oh = 1632 feet is the vertical height y, 

hq =z the tenuity at the end of the interval = 0*956, and therefore the 
average value over the interval is 0*978. 
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The Horizontal Com'ponent. — (4) With the given data the quantity (6) 

becomes 

11-8 

-^ (9) 

76*2 T\p cos ^A 

Assume 1700 an even value of the horizontal velocity. 

The corresponding velocity in the trajectory is 1700/cos 40° = 2219 feet 
per second, for which p is found to be equal to 19. Inserting this value of 
p in (9), together with ti = 0'978, it reduces to O'OIOSG. 

The value of (9) at A, calculated from ^ — 21*03 and t = 0*978 and 
^^ = 40°, is 0*009875. Set these values up at the corresponding horizontal 
velocities (fig. 6), and two points, E and J, are determined on the auxiliary 
curve. 

(5) The area (a) then represents the time for the horizontal velocity to 
change from 1820 to 1700 ; the distance 1700 — M represents this area, and 
M is therefore a point on the time- velocity curve, and OM is the first 
element of the curve. Again the area (b) represents the horizontal distance 
moved through in the time Ot. The length fS is set out to represent this 
distance, so that ISF is a point on the time-distance curve and ON^ is the first 
element of the curve. 

Now set down Oti to represent the time interval T = 1116 second, found 
from fig. 5. Through t\ draw a horizontal to cut the curves, and through the 
points of intersection draw verticals cutting the distance axis in S and the 
velocity axis in 1710. 

(6) Scaling from the curves, OS = 1960 feet, and this is the co-ordinate x ; 
and the velocity at v is 1710 feet per second, and this is the corresponding 
horizontal velocity at B. 

(7) The direction of the tangent at B, 

tan e =.1=: IM^, from which 6 == 39° 45^ 
X 1960 

therefore c^b == 40^-2 (40°~ 39° 45^ - 39|° 

(8) The velocity in the trajectory at B is then 

Y„ — -. oni o = 2201 feet per second, 
^ sm 39i ^ 

1710 

or alternatively V^ = 7~T1FqTo = 2216 feet per second, 

cos oj^ 

a sufficiently near agreement. Point B can now be plotted in fig. 4, and the 
direction of the tangent drawn in. 
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The data at B from which to construct the next element of the trajectory 
are now : — 

Velocity in the trajectory — 2216 feet per second, 
Direction of the velocity ~~3 9 1-° 
Tenuity— 0-978. 

The process is then repeated to get another point C, The second element 
mmi of the time-velocity curve for the vertical component is shown in fig. 5, 
and also the second element mii of the time-height curve corresponding to 
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an assumed vertical velocity at C of 1200 feet per second. The time-interval 
corresponding to this step is 2*146 seconds, and the value of the tenuity is 
found by a vertical from n\ cutting the tenuity curve in qx. Finally, the 
value of the vertical step y\ = 2790 feet. 

The second elements MMi and NIST^ of the time-velocity and the time- 
distance curves for the horizontal component of the motion are shown in 
fig. 6, and the parts of these elements corresponding to the time-interval 
2-146 seconds are X\ = 3440 feet, and the velocity at C is 1508 feet per 
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Upward 



Vertical Components. 



V, 

velocity in 
trajectory. 


V sin (p. 


Angle. 
cp. 


:P' 


T. 

Average. 



ft./sec. 
2375 
2180 
2201 
1888 
1925 
1604 
1676 
1341 
1435 
1076 
1242 

830 
1065 

534 

1000 





11-8 



76'2r^sin<^ + 380 



ft./sec. 


o 


/ 


lbs. 




1525 


40 





21-03 


0-978 


1400 


40 





18 -49 


0-978 


1400 


39 


30 


18-75 


0-897 


1200 


39 


30 


14-65 


0-897 


1200 


38 


36 


15-16 


0-835 


1000 


38 


36 


10-93 


0-835 


1000 


36 


36 


11 -82 


0-753 


800 


36 


36 


7-15 


0-753 


800 


33 


52 


8-76 


0-650 


600 


33 


52 


2-61 


0-650 


600 


28 


52 


5-66 


0-540 


400 


28 


52 


1-08 


0-540 


400 


22 


2 


2-44 


0-450 


200 


22 


2 


0-64 


0-450 


200 


11 


28 


1-86 


0-420 





11 


28 


0-0 


0-420 



-00852 
-00934 
-00987 
-01159 
-01182 
-01450 
-01540 
-01923 
-01900 
-02610 
-02400 
-02940 
0-02870 
-03050 
-03015 
-03100 

Totals . 



Time. 

T. 



Distance, 



sec. 


ft. 


1-116 


1632 


2-146 


2790 


2-630 


2893 


3-463 


8116 


4-510 


3157 


5-340 


2670 


5-920 


1776 


6-115 


611 



31-24 



18,645 



Downward 



Vertical Components. 



y, 




Angle. 




T. 


velocity in 


y sin <p. 




P- 




trajectory. 




cp. 




Average. 



11-8 



76 •2t;^ sin <^- 380' 



Time. 
T. 



Distance. 



ft./sec. 


200 

945 
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1440 
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1190 
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Find 



fi 



ft./sec. 


o / 
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12 I7i 


1-58 


0-48 
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400 


12 m 


14 -54 


0-48 
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0- 


600 
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0- 


600 
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0- 


700 


35 58i 


4-93 


0-804 


0-( 
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nal point 
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T 



-03100 
-03100 
■03210 
•04425 
-03400 
•05685 
•03825 
•05330 
•04200 



Totals 



i 


sec. 


6 


-200 


7 


-635 


9 


-085 


9 


-155 


7 


-340 


39 


•415 



ft. 

620 



2290 



4542 



5950 



5243 



18,645 



second. The co-ordinates of the point C with reference to the point B, fig. 4, 

are therefore 

3/1 = 2790 feet, Xi ^ 3440 feet, 

and the vertical and .horizontal velocities at C are vertical, 1200 feet per 
second, as assumed, and horizontal, 1508 feet per second, as found by the 
process. Finally, the direction of the tangent at C is 
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Flight. 



Horizontal Components. 



Telocity in 
trajectory. 



Y cos ^. 



Angle. 



i?. 



T. 



Average. 



11 '8 



'7Q'2 rp cos ^ 



Time. 
T. 



Distance. 



ft. /sec. 
2376 

2216 

1930 

1686 

1455 

1242 

1090 

1000 



£t./sec. 


o / 
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lbs. 

21-03 

21-03 
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15-24 

11-97 

11-97 

9-12 

9-12 

5-66 

5-66 

2-70 

2-70 

1-86 

1-86 



0-978 
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0-835 
0-835 
0-753 
0-753 
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0-540 
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0-450 
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-00987 
-01075 
-01178 
-01430 
-01578 
-01890 
-02140 
-02630 
-03150 
-04300 
-05780 
-08180 
-13750 
-16400 
0-20250 



Totals 



sec. 
1-116 
2-146 
2-630 
3-463 
4-510 
5-340 
6-920 



6-115 



31-24 



ic. 
1960 
3440 
3686 
4414 

5177 
5607 
5900 
5876 



36,059 



Flight. 



Horizontal Components. 
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Angle. 




T. 


11 '8 


Time. 






velocity in 


T cos (p. 
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JLX O 




Distance. 




76 '2tp cos (() 




trajectory. 




<p. 




Average. 
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ft./sec. 


ft./sec. 


o / 


lbs. 






sec. 


ft. 




942 


942 





1-67 


0-438 


-22500 
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0-438 


-24700 


6-200 


5760 
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916 
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0-48 
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1-58 
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7-340 


5308 




39 -416 


32,958 







tan^ = 



2790 
3440 



0-811, giving e = 39° 3', 



from which 

^c = 39° 30'-2(39° 30'-39° 3') = 38° 36'. 

The details of the whole calculation for the 13 points shown on fig. 8 are 



given in the annexed Table. 



